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ABSTRACT 

The simplest theory describing large-scale redshift-space distortions (RSD), based on 
linear theory and distant galaxies, depends on the growth of cosmological structure, 
suggesting that strong tests of General Relativity can be constructed from galaxy 
surveys. As data sets become larger and the expected constraints more precise, the 
extent to which the RSD follow the simple theory needs to be assessed in order that we 
do not introduce systematic errors into the tests by introducing inaccurate simplifying 
assumptions. We study the impact of the sample geometry, non-linear processes, and 
biases induced by our lack of understanding of the radial galaxy distribution on RSD 
measurements. Using LasDamas simulations of the Sloan Digital Sky Survey II (SDSS- 
II) Luminous Red Galaxy (LRG) data, these effects are shown to be important at the 
level of 20 per cent. Including them, we can accurately model the recovered clustering 
in these mock catalogues on scales 30-200 h~ x Mpc. Applying this analysis to robustly 
measure parameters describing the growth history of the Universe from the SDSS-II 
data, gives f(z = 0.25)ct 8 (z = 0.25) = 0.3512 ± 0.0583 and f(z = 0.37)a a (z = 0.37) = 
0.4602 ± 0.0378 when no prior is imposed on the growth-rate, and the background 
geometry is assumed to follow a ACDM model with the WMAP + SNIa priors. The 
standard WMAP constrained ACDM model with General Relativity predicts f(z = 
0.25)er 8 (z = 0.25) = 0.4260 ± 0.0141 and f(z = 0.37)ct 8 (z = 0.37) = 0.4367 ± 0.0136, 
which is fully consistent with these measurements. 
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1 INTRODUCTION 

The statistical quantification of Redshift-Space Distortions 
(RSD) provides a robust method for measuring the growth 
of structure on very large scales. RSD arise because we in- 
fer galaxy distances from their redshifts using the Hubble 
law: the radial component of the peculiar velocity of indi- 
vidual galaxies will contribute to each redshift and be mis- 
interpreted as being cosmological in origin, thus altering our 
estimate of the distances to them. The measured clustering 
of galaxies will therefore be anisotropic and the additional 
radial signal can be used to determine the characteristic am- 
plitude of the pair- wise distribution of the peculiar velocities 
at a given scale, which in turn depends on the growth rate. 

Many previous analyses have used RSD to measure the 
cosmological growth rate using both th e correlation function 
and power spectrum (see, for example, lHawkins et al.ll2003l ; 
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Percival et al ]|2004l : IZehavi et al.ll2005l;lTegmark et al ]|2006l: 



2008; 



Cabre fc G aztanagal 120091 : ISong et all 



Guzzo et al. 

201 if) . In general these studies used clustering information 



over a small range of scales, and a simplified modelling pro- 
cedure in order to make the measurements. 

Large-scale RSD measurements provide results that 
can be compared to direct measurements of peculiar 
velocities in the local Universe: both observations de- 
pend on the amplitude of the velocity field. Recent 
analyses of the local data seem to indicate the pres- 
ence of unexpectedly large bulk flows, 2a highe r than 
ACDM predictions JWatkins. Feldman fc Hudson) 120091 : 



iFeldman. Watkins fc Hudson! 1201(1 iMacaulav et al] [201ll V. 
although Nusser fc Davis! 1 201lT l present more compatible 
measurements. Large bulk velocities were also detected 
through measurements of the kinet ic Sunyaev-Zeldovich e f- 
fect on the X-ray cluster catalog (Kashlin skv et al] 12009). 
The excess velocities detected are at odds with the previous 
RSD measurements from the 2-degree Field Galaxy Redshift 
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Survey (2dFGRS; IColless et alj|2003h and the Sloan Digital 
Sky Survey (SPSS; I York et alj|2000j > discussed above, which 
give results broadly consistent with ACDM models. There 
is therefore strong motivation for considering if systematic 
effects could be affecting either set of observations. 

If we assume that observed galaxies are sufficiently far 
away that their separations are small compared to the dis- 
tances between them and the observer (the "plane-parallel" 
approximation) then, to linear order, the relationship be- 
tween the redshift-space galaxy power-spectrum P| g , the 
real-space matter power-spectrum P^m an d the growth rate 
is simple (|Kaiserll 19871 ; lHamiltonlll997h . 



tion function and power-spectru m as a function of all three 
variables have been computed (IZaroubi fc Hoffman! Il993l; 
Szalav. Matsubara fc Landvf 19981 ; ISzapu di 2004; Matsubaral 



2004; Papai fc Szapu di 2008). In fact, the wide-angle corre- 



P g s g (fc, M ) = P 1 L m (fc)(6 + / M 2 ) 2 



(1) 



where b accounts for a linear deterministic bias between 
galaxy and matter overdensity fields, / is the logarithmic 
derivative of the growth factor by the scale factor / = 
dlnG/dlna, and fi is the cosine of the angle to the line-of- 
sight. In our paper we study possible theoretical systematics 
beyond the model of Eq. JTJ , that could effect the measure- 
ments of clustering on large scales including the effects due 
to wide-angle corrections, large-scale nonlinearities, sample 
geometry and the effects of the radial model for the distri- 
bution of galaxies. 

Nonlinear effects change the real-space matter power 
spectrum, the velocity power spectrum, the matter-velocity 
cross-correlation, an d introduce furthe r fj, dependent terms 
into this expression l|Scoccimarroll2004h . On small scales the 
dominant nonlinear co ntribution com es from the Fingers- 
of-God (FOG) effect (| Jacksonl Il973 ) . FOG arise because 
within dark matter halos the velocities of galaxies quickly 
become virialized and their power-spectrum is highly non- 
linear. This effect can be approximated by including a 
phenomenologi cal term in Eq. fll tha t reduces power on 
small scales l|Peacock fc Dodda Tl996') or using a more 
complicated expression based on high er order computa- 
tions in perturbation theory (see e.g., IScoccimarrol 1 20041; 
iTaruva. Nishimichi fc Saitol l2010t ). The phenomenological 
dampin g terms used to describe the FOG effects are not ac - 
curate (|Scoccim arro 20041; I Jennings. Baugh fc Pascolill201ll ) 
and the results of perturbation theory are not easy to 
implement in a computationally fast and efficient way. 
The effects of nonlinear growth on the real-space power- 
spectrum are also important and difficult to model for an 
arbitrary cosmological model. Although in principle these 
nonlinear effects can be estimated analytically using per- 
turbation theory, comparison of different perturbation the- 
ory methods to the results of high-resolution N-body sim- 
ulations shows that at low redshifts the range of scales 
where the perturbation theory is reli able is rather small 
l|Carlson. White fc Padmanabhanl2009l ) . In addition, the as- 
sumption that the bias between matter overdensities and 
galaxies is line ar is not accurate even for the scales as large 
as 30/i _1 Mpc l|Okamura fc JindbOllI ). 

Wide-angle corrections are needed because, if the angle 
a that a galaxy pair forms with respect to the observer is 
large, the distance between galaxies is comparable to their 
distance to the observer and the "plane-parallel" approxi- 
mation (and hence Eq. [TJ breaks down. The redshift-space 
correlation function and the power-spectrum in this case will 
also depend on the third variable that could be chosen to 
be the angle a. The wide-angle linear redshift-space correla- 



lation function does not deviate significantly from its "plane- 
parallel" counterpart if the opening angle a is less than 10° . 
In previous wor k we validated this work by analysing m ock 
galaxy catalogs (Raccanc lli. Samushia fc Perciva3l2010t ). 

For surveys that cover a significant fraction of the sky, 
the distribution of galaxies pairs becomes non-trivial. The 
survey geometry results in the galaxy pair distribution that 
has a complicated dependence on the variables r, /j, and a, 
since not all sets of their combinations are equally likely 
or even geometrically possible. In particular the distribu- 
tion of fi does not correspond to that of an isotropic pair 
distribution. This will strongly bias the measurement of 
angular momenta of the correlation function, and in fact, 
often dominates over differences between "plane-parallel" 
and "wide-angle" effects for galaxy pa irs with the same {i 
l|Raccanelli. Samushia fc Percivalll2010l ). 

RSD data on very large scales, although in principle 
available in current data sets, do not contribute significantly 
to current data analyses. The reason for this is twofold: the 
signal-to-noise of currently available clustering data becomes 
low at scales larger than 100 h~ x Mpc so most of the avail- 
able cosmological information is on smaller scales; also the 
large scale clustering measurements are vulnerable to dif- 
ferent observational (improper modelling of seeing, galactic 
extinction, etc.) and theoretical systematic effects which, if 
not taken into account properly, could strongly bias results 
of data analysis. In particular, our ability to model the ra- 
dial galaxy distribution accurately can cause strong effects 
on these large-scales, and is worthy of further investigation 
l|Percival et al.ll2010l ; lKazin et alJb'ulCh . Being able to model 
these data has many advantages. First, if accurate measure- 
ments are available, more data will result in stronger con- 
straints on cosmological parameters. In addition, measure- 
ments on large scales are significantly less affected by the 
systematics introduced by nonlinear phenomena. Some im- 
portant physical processes can be measured only on very 
large scales. For example, non-Gaussian initial conditions, 
if present, will affect the rea l -space galaxy clustering o n 
large scales l|Dalal et al l 120081 : iDesiacques fc Seijakl |201Ch . 
and could be compared against the RSD signal, which de- 
pends on the matter field. 

We investigate the significance of these effects by per- 
forming an analysis on a large suite of N-body simulations, 
testing for systematic effects that could result in real data 
giving a signal different from the plane-parallel linear RSD 
formula. Using mock samples, we are able to accurately 
fit the expected correlation function on scales between 30- 
200/i _1 Mpc, to a level well below the statistical error on 
the measurement from any one sample. 

We apply the knowledge learned in this analysis to ro- 
bustly measure RSD in the Sloan Digital Sky Survey (SDSS) 
Data Release 7 (DR7) sample of Luminous Red Galax- 
ies (LRGs) and measure cosmological parameters describ- 
ing amplitude and growth of the perturbations in differ- 
ent models. We find that the accuracy of SDSS DR7 data 
is at the threshold where the inclusion of RSD informa- 
tion on scales larger than 100 ft -1 Mpc affects the measure- 
ments but does not improve the result significantly. In ad- 
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Figure 1. Histogram showing the redshift distribution of galaxies 
in the SDSS DR7 LRG catalog used in our analysis. 



dition, we find that the non-linear and survey geometry ef- 
fects are significant for this sample. The next generation of 
spectroscopic surveys, such as Baryon Oscillation Spectro - 
scopic Sur vey (BOSS) JSchlegel, White fc Eisenstem[|2009l) 
BiqB OSS (|Schlegel et all 120091 ) and Euclid l|Laureiis etall 
2009) will enable us to measure the clustering of galaxies at 
large scales with extremely high accuracy and using RSD 
data on large scales will be crucial for constraining growth 
of structure in the Universe and the nature of gravity in a 
robust way. On large scales, ignoring effects such as wide- 
angles would lead to systematic deviations. However, these 
effects are known, so can be easily included. On small scales, 
non-linear behavior is harder to model, so there may be un- 
avoidable additional errors here. 

This paper is organised as follows, in Sec. [5] we describe 
our data set and the method of measurement. In Sec. [3] we 
present our theoretical model for the RSD on very large 
scales and check its validity with the results of N-body sim- 
ulations. In Sec. [4] we describe theoretical models and as- 
sumptions made when extracting cosmologically relevant in- 
formation. In Sec. [S]we present results of our data analysis. 
In Sec. [6] we conclude and discuss the relevance of the con- 
tents of our paper to future spectroscopic galaxy surveys. 



2 CALCULATING MOMENTA OF THE 
CORRELATION FUNCTION 

2.1 The SDSS data 

We use data from the SDSS data rel ease 7 (DR7), wh ich 
obtained wide-field CCD phot ometry (|Gunn et al]|l998l ) in 
five passbands (u, g, r, i, z\ e.g. jFukugita et al.lll996h . amass- 
ing nearly 10,000 square degrees of imagin g data for which 
objec t detection is reliable to r ~ 22 (|Abazaiian et al] 
2009). From these photo metric data, Luminous Red Galax- 
ies (LRG) were targeted |Eisenstein et al]|200ll ) and spec- 
troscopically observed, yielding a sample of 106,341 LRGs 



in the redshift bin 0.16 < z < 0.44. The redshift distribution 
of galaxies in this catalog is shown on Fig. [T] 

To study clustering properties of LRGs we create a ran- 
dom catalog that has unclustered "galaxies" randomly dis- 
tributed with the same angular mask as SDSS DR7. The an- 
gu lar distribu ti on of these galaxies was chosen as described 
m iReid et all (|2009h . The method for, and effect of esti- 
mating the expected radial distribution of the galaxies is 
described in Section 12.31 Our random catalog has approxi- 
mately 50 times more objects than the real catalog. 



2.2 Methodology 



We assign each galaxy a weight 
1 



1 + n(r)P 



(2) 



where n(r) is a local density of galaxies in units of 
(ft -1 Mpc) in a neighbourhood of the galaxy of interest 

located at a position r and P = 10000 ( h^ 1 Mpc) 3 . 
This weighting is optimal for the premise that 
;alaxies Poisson sample t h e un derlying matter field 
Feldman. Kaiser fc Peacockl Il994h . Recent work has 
shown that it may be possible to beat this if we 
can estimate the mass associate d with each galaxy 
l|Seliak. Hamaus fc Desiacauesl 120091 ') : we do not attempt 
this in our work. 

In order to extract information about the evolution of 
structure growth, we divide the LRG sample in two redshift 
bins so that the weighted number of galaxies 



AT. 



(3) 



is approximately equally split between them. The two red- 
shift bins are 0.16 < z < 0.32 and 0.32 < z < 0.44 while the 
effective redshifts for our two bins are z — 0.25 and z = 0.37. 
Effective redshift is defined as 



ZcB 



N W (N W - 1) 



y^y^«;jWj(zj+zj), 



(4) 



where N w is given by Eq. Q. 

For each pair of objects in our galaxy catalog, random 
catalog or cross pairs between galaxy and random catalogs, 
we compute the distance between objects r, the angle a 
that the objects make with respect to us, and the cosine 
of the angle that the bisector of the angle between the ob- 
jects makes with the line connecting them (assuming a flat 
geometrical model) /x. Due to statistical isotropy about the 
observer, these three variables are sufficient to completely 
describe the RSD expected for each pair. We bin r in 65 
equal logarithmic bins from 1 Mpc to 200 ft -1 Mpc, \i in 
200 equal bins from to 1, cos(a) in 400 bins from -1 to 1 
and count the number of galaxies in each bin0 



1 Later in the paper we will show that for this particular observed 
geometry wide-angle effects are negligible, which means that for 
this particular case the a-label could have been drop from the 
beginning. We still keep the a-label in the rest of this paper for 
the completeness of formalism. 
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To convert angular and redshift separations of galax- 
ies into physical separations a fiducial cosmological model 
is needed. We compute distances in a spatially-flat ACDM 
fiducial model. If the real geometry of the Universe is dif- 
ferent from the one described by our fiducial model it will 
bias the measurements of clustering through the Alcock- 
Paczynski effect. We discuss this issue in Sec. 13.51 

In the plane-parallel approximation, all of the available 
linear RSD information can be extracted from the zeroth, 
second and fourth Legendre momenta of th e galaxy correla - 
tion function with respect to the variable /i l|Hamiltonll 19921 ). 
Given that we expect that wide-angle and non-linear effects 
will give relatively small deviations about this approxima- 
tion, we should expect that, even in the more general case, 
these momenta contain almost all of the available RSD in- 
formation. We therefore choose to fit to these measurements 
in our work. This will be discussed further in Section [3] To 
estimate those three w e use Landy-Szalay type estimators 
l|Landv fc Szalavlll993h 

6( r = ^ {[DD(r l ,fi j ,a k ) - 2DR(n, fj,j, a k ) 

+ RR(n, fij , a k )] P t ) } / ^ RR(n, ^ , a k ) , (5) 

i,k 

where DD(n, fij, a k), DR(n, /ij , a fe ) and RR(ri, fij, are 
the numbers of galaxy-galaxy, galaxy-random and random- 
random pairs in bins centered on Vi, fij and at- The Pi are 
the I th Legendre momenta. 

The measurements of £e(r) from all LRGs in our cat- 
alog, together with theoretical predictions of spatially-flat 
ACDM model with fi m = 0.25 are shown on Fig. Note 
that, Fig. [2] shows only statistical errors and does not show 
extra systematic errors due to uncertainty in radial distribu- 
tion of galaxies (for details on systematic errors see Sec. 12. 31) . 
Also the statistical errors are computed from diagonal ele- 
ments of the covariance matrix only, the whole structure of 
the covariance matrix is such that measurements are more 
likely to be systematically above or below theoretical line 
rather then randomly scattered around the theoretical pre- 
diction as for noncorrelated Gaussian variables. 

RSD measurements are also often extra cted from the 
normalised quadrupole Q(r) (|Hamiltonll 19921 ) . defined as 

Q(r) = ^ ■ (6) 

£o(r) - ^ / £o(r>' 2 dr' 
Jo 

We can form an estimator for Q(r) by replacing integral in 
Eq. |6) by a discreet sum 

$('«) = 1^ . (7) 

3=0 

The measured Q(r) from the SDSS DR7 LRG data is shown 
in Fig. [3] The details of how the statistical error bars are 
computed are discussed in Sec. 13.71 
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Figure 2. Measurements of £ e (r) from SDSS DR7 LRGs in a 
redshift range 0.16 < z < 0.44. The statistical error-bars were 
calculated as described in Section \3\ and represent only the diago- 
nal elements of the whole covariance matrix. The absence of lower 
error-bar on some measurements indicates that they are consis- 
tent with zero. Solid line shows a theoretical prediction with the 
shape corresponding to the best fit cosmology to current WMAP 
and SNIa measurements and the amplitude given by the best-fit 
values to the data. 
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Figure 3. Measurements of Q(r) from SDSS DR7 LRGs in a 
redshift range 0.16 < z < 0.44. Error-bars were calculated as 
described in Section [3] and show only the diagonal elements of 
statistical covariance matrix. Solid line shows a theoretical pre- 
diction with the shape corresponding to the best fit cosmology to 
current WMAP and SNIa measurements and the amplitude given 
by the best-fit values to the data. 



2.3 Modelling the redshift distribution of SDSS 
LRGs 

To measure a correlation function from a survey accurately, 
we must know what an unclustered distribution of galaxies 
would look like in the same volume. The unclustered distri- 
bution can in principle be derived by averaging observations 
in many unconnected regions. Since the real data covers only 
a relatively small volume the expected galaxy density (in the 
absence of clustering) is hard to determine in this way (cos- 
mic variance). The wrong estimate of unclustered distribu- 
tion will bias the measurements of the correlation function. 
This effect is especially important on large scales where the 
fluctuations we wish to measure are small. 

In this paper we compute correlation function by using 
a spline fit to th e galaxy redshift d istribution (with param- 
eters as given in lPercival et al.M2010r ) . We will refer to it as 
a "spline" random catalog. The exact form of the random 
catalog will depend on the position and number of nodes 
used for the spline. In the limiting case when the number 
of nodes goes to infinity while the spacing between nodes 
goes to zero we will have a random catalog that has ex- 
actly the same redshift distribution of galaxies as data: in 
effect we assign a randomly chosen galaxy redshift to the 
random object . We refer to it as a "z-shuffled" catalog. 
We also construct a random catalog by randomly mixing 
angular positions and redshifts in the galaxy catalog (later 
referred to as "3D-shuffled" catalog). We should expect the 
shuffled catalogs to remove some structure, as fluctuations 
in the galaxy density caused by large-scale structure will 
be smoothed. The "spline" and "z-shuffled" catalogs, unlike 
the "3D-shuffled" catalog, have angular positions of objects 
choosen at random within the sample angular mask. 
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Figure 4. Radial distribution of galaxies in different random cat- 
alogs. Dashed line corresponds to the "shuffled" random catalog, 
which has radial distribution identical to that of data. Dashed line 
corresponds to "spline" random catalog and solid line corresponds 
to "proper" random catalog. 



To quantify the possible systematic offset induced by 
improper modelling of the radial distribution of galaxies we 
use large suite of LasDamas N-body simulations (McBride 
et al., in prep) which are designed to replicate the observed 
geometry of the SDSS-II (for more details on LasDamas sim- 
ulations and how we use them see Sec. I3.7[> . For the mock 
catalogs the unclustered redshift distribution of galaxies is 
known and we will refer to the random catalog based on 
this known distribution as a "proper" random catalog. We 
compute correlation function of mocks using the "proper", 
the "spline" , the "z-shuffled" and the "3D-shuffled" random 
catalogs. The radial distribution of galaxies in each of these 
random catalogs, for one of the LasDamas mock, is shown on 
Fig-S] "z-shuffled" and "3D-shuffled" catalogs have identical 
redshift distribution of points but different angular distribu- 
tion. 

Figure [4] shows that although the "spline" random cat- 
alog follows the general shape of the "proper" random cat- 
alog, it does not reproduce the real radial distribution of 
galaxies accurately. The induced systematic errors on mea- 
surements of £o and £2 when using different random catalogs 
catalog with respect to the "proper" catalog, averaged over 
all 80 LasDamas mock, are shown on Fig. [5] along with the 
statistical errors. 

Figure [5] shows that when the "3D-shuffled" catalog is 
used the systematic offset is larger then statistical errors on 
all scales. The statistical errors are larger for £2 compared to 
£o- This is not surprising since the errors in redshift distri- 
bution affect clustering mainly in radial direction. £0 mea- 
sures average clustering in all direction, while £2 measures 
the excess of clustering across the line-of-sight compared to 
radial and is expected to be affected by this systematics 
more. For the "z-shuffled" and "splined" random catalogs 
the systematic offset is a fraction then statistical errorbars 
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Figure 5. Relative errors introduced in the measurement of the 
first two even Legendre momenta of the correlation function using 
different random catalog. The solid line shows statistical errors. 



and "z-shuffled" catalog seems to be performing better com- 
pared to the "splined" catalog. 

In our analysis we use the "splined" catalog to compute 
correlation functions and ignore this systematic offset since 
it is small compared to current errorbars. The exact scal- 
ing of the offset with galaxy number density, and radial fit 
technique will be investigated in a separate paper. 



2.4 Excess of power on large scales 

The top panel of Fig. [2] shows that there is an excess of 
power in the measured monopole of the correlation function 
with respect to the predictions of flat WMAP 5 normalised 
ACDM model. This excess has been observed previously 
in spectroscopic (|Okumura et al . 2008: ISanchez et al.ll2009l : 



ICabre fc Gaztanagal I2009I; iKazin et all |2010l) and photo- 
metric (iThomas. Abdalla fc Lahavll201ll : ICrocce et aT]l201ll : 
iRoss et al.ll201ll ) data sets. Recent analysis of mock catalogs 
resembling SDSS DR7 showed that 6 out of 160 measured 
correlation functions wer e positive up to t he scales of 200 
/t^Mpc (see, Fig. 12 in IKazin et~aTl l201fj| 'l . This suggests 
that the excess seen is only mildly statistically significant. 

If the signal is physical, the modifications to the stan- 
dard cosmological model that result in stronger clustering, 
such as some modifie d gravity theories or clustering dark 
energy (Takada 200(|), could explain this anomaly. Other 
possible explanati ons are the presen ce of large non-Gaussian 
initial conditions l|Dalal et al.ll2008h or isocurvature pertur- 
bations. It should be noted however that the amount of non- 
Gaussianity required to generate such a big offset is ruled out 
by ot her observations (see, e.g., iDesiaccmes. Seliak fc Ilievl 
2010). Possible observational systematics include, for in- 
stance, improper modelling of extinction and seeing. Both 
of these could introduce spurious extra angular fluctuations 
in the da ta that would l ater b e misinterpreted as an excess 
of power. ISanchez et all (|2009h showed that if there is a sys- 
tematic constant shift in measured correlation function (due 
to calibration errors or evolutionary effects) this does not 
bias the estimated best-fit valus of cosmological parameters 
significantly. 

Inaccuracies in the measured clustering induced by as- 
suming an incorrect radial distribution, are also large an d 
could in principle explain the anomaly |Kazin et alj|20ich . 
The error rescaling suggested in previous section makes the 
inconsistency with the standard model on very large scales 
less severe □ 



3 MODELING RSD ON LARGE SCALES 
3.1 Plane-Parallel, Linear model 

The linear, pl ane-parallel m odel for RSD is often termed the 
Kaiser model (jKaiserll 19871 ). In the following, we follow stan- 
dard convention and denote the observed galaxy overdensity 
field by S s , with a superscript s in redshift-space and r in 
real-space. A given Fourier k mode of this overdensity can 
be expressed to linear order in overdensity as 



S s e (k) = S r g (k)-^ 2 9 s (k), 



(8) 



where p is the cosine of an angle with respect to the line of 
sight and 9 g = V • u is a divergence of the galaxy velocity 
field. We follow the commonly adopted assumption that this 
is equal to the divergence of the matter field, assuming no 
velocity bias, i.e. we assume that b v (k) — 1, where # g (fc) = 
b v (k)0m(k). A subscript g shows that a quantity relates to 
the galaxy field, and a subscript m denotes the matter field. 

The two-point function of this overdensity field is 
anisotropic, 



2^P r s g{k)+^P r e , 



(k), 



(9) 



2 Also note that the measurements of spherically averaged cor- 
relation function at different scales are strongly correlated. This 
makes the deviations to the one side of the model prediction more 
probable, the significance of this deviation being smaller then 
what it would be for uncorrelatcd measurements. 
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where P xy = (5 x S y ) denotes a cross-power-spectrum of fields 
x and y. 

If we further assume that, in real-space, the overdensi- 
ties in the galaxy field are linear functions of overdensities in 
the matter field <5 g = 6<5 m and the velocity divergence can be 
related to the matter overdensities using the linearized con- 
tinuity equation Sg — —f8 m , then the redshift-space power- 
spectrum can be simply expressed in terms of the real-space 
power-spectrum as 



p g s g (fc,/i) = (b+fi/yp^ik). 



(10) 



The proportionality constant b between matter and galaxy 
overdensities is the bias factor and the coefficient / between 
velocity divergence and the matter overdensity is equal to 
the logarithmic derivative of the growth factor by the scale 
factor dlnG/dlna, which follows from the continuity equa- 
tion combined with scale-independent growth. 

The redshift-space correlation function is given by the 
Fourier transform of Eq. <(9j 



Pg g (fc,/i) exp(— ikr)d k, 



(11) 



and can also be expres sed in terms of its Legendre momenta 
l|HamiltorJll992l . ll997f ). 



3.2 Wide-angle effects 

When the distance between galaxy pairs is compara- 
ble to the distance between galaxies and the observer, 
the theory of Sec. 13.11 can not be used to describe 
RSD effects. The redshift-space correlation function be- 
comes a function of three variables that can be chosen 
to be the separation between galaxies r and the two 
angles cj>i and <f>2 that galaxies form with an arbitrary 
z ax is with respect to the observ er (Zaroubi & Ho ffman] 
1993!: ISzalav. Matsubara fc Landvl Il998l : ISzapudil |2004 ). 
Papai fc Szapudil (120081 ) showed that this correlation func- 
tion, when expanded in tripolar spherical harmonics, 



£s(r,(f>i,ct>2) = 



E 



(r ) </>i,^2)a i ^ 2 ,f(xi,x 2l x), (12) 



gives only a few non-zero terms in the absence of an obser- 
vational window. Here xi and X2 are the unit vectors in the 
direction of two galaxies and x is a unit vector pointing in 
the direction from galaxy one to galaxy two. 

Eq. (|12|l can be recast as a function of variables r, fj, and 
a, where a is an angle the galaxies make with respect to the 
observer. This set of coordinates is invariant with respect 
to rotation and more straightforward to use in data analy- 
sis. There are three reasons for differences between "plane- 
parallel" and "wide-angle" predictions: 

(a) The "wide-angle" correlation function £(r, /x, a) de- 
pends on a, while the "plane-parallel" one doesn't; 

(b) The coefficients B ll ' i2,1 depend on the density of 
galaxies n(z) as a function of redshift. This implies that 
the RSD effect will depend on the spatial distribution of 
observed galaxies; 

(c) The distribution of galaxies in fi will be non-trivial, 
with some values of [s, not permitted for non-zero a. As a 
consequence, we will not be able to measure pure Legendre 




a (rad) 

Figure 6. Normalised distribution of pairs in SDSS DR7 LRG 
catalog as a function of fi and a at different scales. Top panel 
shows distribution of fi for a fixed r summed over all values of a. 
Bottom panel shows distribution of a for fixed r summed over all 
values of /i. 



momenta of the correlation function, but instead will use 
weighted integrals and biased momenta. 



In iRaccanelli. Samushia fc Percival (|2010h we used simula- 
tions to demonstrate these effects, showing that they should 
be carefully taken into consideration in order to fit the mea- 
sured wide-angle correlation function. In the following we 
describe (a) and (b) as wide-angle effects, whereas (c) is 
termed the "/i-distribution" as it could be applied to plane- 
parallel and wide-angle theory of individual line-of-sight. 

Allowing for full distribution of galaxy pairs, the esti- 
mates of Legendre momenta given by Eq. (J5]) correspond 
to 
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£(r, (j,, a)W(r, fi, a)Pi(p)dfida, 



(13) 



where £(r, a) is given by either the wide-angle formula 
in Eq. (|12|l or its plane-parallel equivalent computed from 
Eq. pop . W(r, n, a) is a weight factor that gives the relative 
number of pairs in a survey that form angles fj, and a for a 
given scale r. The weight factor W is normalised so that 



W(r, (i, a)d/j,da = 1 



(14) 



for all scales r. Ignoring (a) corresponds to setting 
W(r,[i,a 0) = and ignoring (c) corresponds to setting 
W(r,fj,, a) = 1. 

In practice, W(r, fi, a) weights can be computed from 
the random catalog; they will be given by properly nor- 
malised RR(r,(i,a) number counts. In Sec. 12.31 we showed 
that uncertainties in radial distribution (that translate into 
uncertainties in RR counts and consequentially into uncer- 
tainties in the W weights) significantly affect measurements 
of correlation function. This uncertainty will also bias our 
modelling of /^-distribution effects. This effect is of higher 
order and we do not investigate it further in our paper. 

Fig. [6] shows the normalized distribution of pairs in fi, 
and a for different scales for the SDSS DR7 LRG catalog. 
When the a distribution tends towards a delta function cen- 
tred at a = 0, the wide-angle effects (a) becomes negligi- 
ble. When the distribution in fj, tends towards a uniform 
one, (c) becomes negligible. In general, the relative impor- 
tance of the wide-angle effects on the measured correlation 
function depends on the geometry of the survey, its redshift 
range and what scales are considered. The effect is stronger 
for lower redshifts and becomes increasingly important on 
larger scales. Top panel of Fig. [S] implies that for SDSS DR7 
observed geometry it is easier to fit galaxy pairs across and 
along the line of sight rather than for angles in between, 
more so for larger scales. In Sec. 13.71 we will show that for 
the SDSS DR7 geometry, the difference due to (a) and (b) 
in the list above are much smaller than statistical errors and 
can be safely ignored even for scales as large as 200 Mpc. 
The differences due to a nontrivial /^-distribution (item (c) 
in the list above) are larger than wide-angle effects but small 
compared to current statistical errors. They are, however of 
order of few percent at larger scales and will be important 
for future surveys. 



3.3 Nonlinear effects 

The following nonlinear effects, if they are comparable to 
the measurement errors, can make Eq. (|1UJ1 unsuitable for 
analysing RSD data 

(a) Nonlinear contributions to the relationship between 
matter and galaxy overdensities 8^ — h(5^), where h is an 
arbitrary function. 

(b) The relationship between the velocity divergence and 
matter overdensities # g = —f5 m relies on scale-independent 
linear growth coupled with the continuity equation. Also the 
galaxy velocity divergence field must be an unbiased tracer 
of the matter velocity divergence field, i.e., 8 e — # m . This 
formula will break down if these conditions are not met. 



(c) The matter power spectrum itself goes non-linear, be- 
cause of the scale-dependent non-linear growth on smaller 
scales. 

(d) The real-space to redshift-space mapping includes 
higher order terms involving 5g and # g |f| 

In the following, we are only interested in the signal on large- 
scales where linear theory should be strongest. We therefore 
assume that all non-linear effects are small except for (c), 
where we allow the overall powe r spectrum sha p e to d eviate 
from the linear form (but see, iReid fc White] HoTJ)). We 
use mock catalogues in Section ^. 7l to confirm the validity of 
this assumption. 

In order to approximate the non-linear power spectrum, 
we adopt a two-component model, which splits P(k) into a 
"smooth" part that describes the overall shape and a "wig- 
gled" part that describes the Baryon Acoustic Oscillations 
(BAO), 



PbB.o(k, fx) = Pfuii(fc, n) 



(15) 



The "smooth" part is defined by taking some reasonably 
spaced points ki and then interpolating the linear power 
spectrum values betw een those nodes u sing a bi-cubic spline 
interpolation routine (|Press et al.ll 19921). In this work we use 



ki spacing similar to iPercival et al.l |2010| ; we place nodes 
at k = 0.001, k = 0.25 and k = 0.25 + 0.05n where n is 
large enough for the purposes of recovering the correlation 
function by the means of a Fourier transform 

The primary non-linear effect on the BAO compo- 
nent of the power spectrum is a damping on small scales, 
whic h can be well approximated by a Gaussian smooth- 
ing dBharadwail Il996l; ICrocce fc Scoccimarrd l200d 120081 ; 



lEisenstein, Seo fc Whitell2007t lMatsubarall2008al lbl) 



Pbno(k, M) — Pba,o{k;jj) X 



exp 



(1-M 2 )S± | M 2 S|, 



(16) 



where Ex = EoG and En = EqG(1 + /)• So is a constant 
phenomenologically describing the nonlinear diffusion of the 
BAO peak due to nonlinear evolution. From N-body simula- 
tions its numerical value is of order 10 ft -1 Mpc and seems to 
depend linearly on as but only weakly on k and cosmological 
parameters. 

Next order non-linear effect results in a tilt of correla- 
tion functio n on large scales just before the BAO peak (For 
details see, ISanchez et al.l [20091 ). We do not consider this 
and other higher order terms in our computations. Robust 
data analysis of future high quality measurements should 
also include a modelling of this small scale nonlinear effects. 

3.4 Fingers of god effect 

Within dark matter haloes the peculiar velocities of galax- 
ies are highly non-linear. These velocities can induce RSD 

3 This refers to the nonlinear effects in the large-scale coherent 
motions and not to FOG nonlincarities due to virialization of 
galaxy velocities within halos which will be discussed in Sec. 13.41 

4 When making actual Fourier transform we are using much 
denser set of nodes compared to the ones used to define "smooth" 
power-spectrum to make sure that the noise introduced by dis- 
creetness of Fourier transform is small. 
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that are larger than the real-space distance between galax- 
ies within the halo. This gives rise to the observed fingers of 
god (FOG) e ffect - strong e longation of structures along the 
line of sig ht (|jacksonlll973 ). The FOG effect gives a sharp 
reduction of the power spectrum on small scales compared 
to the predictions of the linear model, and is usually mod- 
eled by multiplying the linear power-spectrum by a function 
F(a v ,k, fi), where a v is the average velocity dispersion of 
galaxies within the relevant haloes. The function F is cho- 
sen so that it is small on small scales and approaches unity 
for scales larger t han l/a v . The two most frequently used 
functions are (e.g. lCole et~ai1l 19951 ; [Peacock fc Doddsj|l996h 



fcf, 



M = 



-FLorcntzian(fc,/i 2 ) = [l + (fctf^) 2 ] 
^Gaussian^,^) = exp \ - (ka v /J,) 2 ] . 



(17) 
(18) 



Note that this model is constructed by a rather ad-hoc splic- 
ing of the FOG signal together with the linear model and 
ignores the scale-dependence of the mapping between real 
and redshift-space separations (jFisherl ll995l ; IScoccimarrol 
|2004| . and references therein). In addition, the exact form 
of F(k,/j, 2 ), and the v alue of a v is strongly dependent o n 
the galaxy population ijjing fc Borner 2004; iLi et aLll2007h . 

The Gaussian smoothing in Eq. I|16[l. amongst other 
nonlinear effects, also accounts for the damping due to ran- 
dom velocities described by Eqs. (|17p - (|18|) . In our analysis 
we will use a model given by Eq. (|16[1 that partially includes 
the FOG effect on large scales and will ignore FOG effects 
on small scales. 



3.5 Degeneracy with Alcock-Paczynski effect 

The positions of galaxies in our catalog are given in terms 
of the angular positions and redshifts. To convert angu- 
lar and redshift separations into physical distances the an- 
gular and radial distances as functions of redshift are re- 
quired. Those functions depend on the adopted cosmolog- 
ical model. We perform our pair count assuming a fidu- 
cial, spatially-flat cosmology with fi m = 0.25. If the real 
cosmology is significantly different from the fiducial one, 
this difference will introduce additional anisotropies in the 
corre lation function through th e Alcock-Paczynski (AP) ef- 
fect l|Alcock fc Paczvnskilfl979l'). This can significantly bias 
the measurements of growth dBallinger. Peacock fc Heavens! 
ll99rj ; ISimpson fc Peacockll2010l ; ISamushia et alj|2010h . 

In the presence of the AP effect the redshift-space 
power-spectrum is 



P B (k'.. 



fi ,QX,Q||,p) 




12 r 
V / 



F 2 + m 2 (1- F 2 ) 



1 + ^ 2 (^-l 



(19) 



where p are standard cosmological parameters determining 
the shape of the real-space power-spectrum, k' and // are the 
observed wavevector and angle, related to the real quantities 
by 



Ml 



ot\\k\\ 
a±_k± 



(20) 
(21) 



k' 2 + k' 2 



(22) 



the a 1 1 and a± are the ratios of angular and radial distances 
between fiducial and real cosmologies 



H 



fid 



(23) 
(24) 



and F = a\\/a±. 

Ignoring the AP effect is equivalent to assuming that a 
factors are equal to unity in Eq. (|19|l . This assumption can 
bias estimates of growth parameters and their uncertainties. 

We estimate to magnitude of this effect for our analysis 
using Fisher matrix method. We c ompute a Fisher matri x 
for the SDSS-II like survey following ISamushia et al.l (|2010h . 
This Fisher matrix is an optimistic estimate of the inverse 
covariance matrix on the parameters 6, /, an, a± and p. 
The covariance matrix is an inverse of this fisher matrix. 
Ignoring the Alcock-Paczynski effect is equivalent to remov- 
ing rows and columns corresponding to an and a± first, as 
if they were perfectly known, and only then inverting the 
fisher matrix to get covariance of b and /. The more accu- 
rate approach is to invert the Fisher matrix directly without 
assuming that the as are known. 

In our data analysis, we will apply a prior based on the 
WMAP and SNIa data on the background geometry of the 
Universe (see Sec. HJ. To reflect this in our Fisher matrix 
computations we first add this prior to the Fisher matrix 
elements corresponding to an and ct±_ and only then invert 
the whole matrix to get covariances on b and /. We compare 
the result with the resulting covariances when the AP effect 
is ignored. 

To do this we use MCMC chains corresponding to 
WMAP and SNIa joint constraints on spatially-flat WCDM 
Universe from WMAP LAMBDA websit<0 to estimate 
Fisher matrix of parameters fi m , wo and Ho . The actual pri- 
ors are fi m = 0.276 ± 0.020, w = -0.969 ± 0.054, w a = 0, 
fik — while the errorbars are Gaussian and slightly cor- 
related. We transform this into a Fisher matrix on an and 
a±. 

Figures [7] and |S] show the effects of AP on the measure- 
ments of different growth parameters (for the description 
of parameters 7 and / see Sec. [4}. These figures show that 
in general ignoring the AP effect results in gross underes- 
timation of the error bars. Real uncertainties on 7 and / 
are few times larger then what we would get when ignoring 
AP. After applying the strong prior on the background ex- 
pansion, however, almost all of this degeneracy is removed 
and the uncertainties in the measurements of growth and 
bias are almost iden tical to the case with no AP, consis- 
tent with the work of lSamushia et al] |2010), which showed 
the importance of model assumptions on this measurement. 
We conclude that, for the models we test, the effects of the 
degeneracy between RSD and AP on the error bars of our 
measurements are very small and can be safely ignored pro- 
vided we adopt joint WMAP and SNIa priors. 



http : //lambda . gsf c . nasa . gov/product/map/dr4/params/wcdm_sz_lens_wmap7_ 
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Figure 7. Fisher matrix predictions of correlated constraints on 
parameters b and 7 from a SDSS-II like survey. The dashed con- 
tours correspond to the most general case when the AP effect 
is not ignored, the dotted lines correspond to the case when the 
AP effect is ignored and the solid lines correspond to the case 
when the AP effect is ignored but a strong prior is put on the 
background cosmology. The solid and dotted lines are almost in- 
distinguishable by eye on this plot. 
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Figure 8. Fisher matrix predictions of correlated constraints on 
parameters b and / from a SDSS-ff like survey. The dashed con- 
tours correspond to the most general case when the AP effect 
is not ignored, the dotted lines correspond to the case when the 
AP effect is ignored and the solid lines correspond to the case 
when the AP effect is ignored but a strong prior is put on the 
background cosmology. The solid and dotted lines are almost in- 
distinguishable by eye on this plot. 
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Figure 9. Relative impact of nonlinear and wide-angle effects 
compared to the statistical errors on the measurements of !;e(r) 
from SDSS DR7 data. 

3.6 Relative importance of different effects 

We now consider the relative importance of the modifica- 
tions to the linear plane-parallel model, described in pre- 
vious subsections, as a function of scale. We have assumed 
a spatially-flat ACDM model with f2 m = 0.25 throughout. 
Figs.|9]and[l0]compare statistical errors on measurements of 
t;i(r) and Q(r) (for the details of how these statistical errors 
are estimated see Sec. I3.7fl with the differences in best-fit 
theoretical models calculated with or without the modifica- 
tions considered above. 

To estimate the impact of different systematics we first 
compute a theoretical correlation function for our fiducial 
model £(r) u11 including all effects. We compute linear £(r) 
using CAMB l|Lewis, Challinor fc Lasenbvll200b1 ). Then we 
recompute the same correlation function by ignoring each 
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Figure 10. Relative impact of nonlinear and wide-angle effects 
compared to the statistical errors on the measurements of Q(r) 
from SDSS DR7 data. 



of the systematic effects in turn to see by how much this 
changes our theoretical estimates at different scales. Nonlin- 
ear diffusion of the BAO peak is modelled with the Eq. ()16[) 
with E = 8 ft" 1 Mpc0 the FOG effect with the Gaussian 
damping function of Eq. (|18[) with a v — 3.5 h" 1 MpcQ the 
effect of jj, distribution is studied first by using the real 
distribution of angles in SDSS geometry for W and then 
assuming it to have a uniform probability over all angles; 
the magnitude of the effect to small scale nonlinearities 
is calculated by comparing correlation functions computed 
from linear power-spectrum to the nonli near one (compute d 
using the HALOFIT fitting formula of ISmith et aHl2003tl ; 
the wide-angle effects are estimated by substituting the full 
wide-angle correlation function by a two-dimensional plane- 
parallel approximation. 

The results are shown on Figs. [9] and \W\ The FOG ef- 
fect and corrections to the shape of the correlation function 
due to nonlinear growth of structure are only important on 
smaller scales and are lower than the measurement errors 
on the scales larger than 20 ft -1 Mpc for the monopole and 
scales larger than 30 h" 1 Mpc for the quadrupole. It should 
be noted that the FOG induced relative errors as shown on 
Figs. [5] and [TD] correspond to the difference between using 
the anisotropic Gaussian dumping of Eq. (|18[) and ignoring 
the effect altogether, for a wrong but reasonable FOG mod- 
elling the induced systematic errors will be smaller. The wide 
angle effects appear at the scales of about 70 h~ l Mpc, as an- 
ticipated based on Fig. [6] but are less then 1 percent even 



6 The value is consistent with the estimate of So for a real- 
space power-spectrum on low redshifts and at large scale 
computed from N-body sim ulations in a standard ACDM 
llEisenstein. Seo fc Whitell2007r> . 

7 This val u e is consistent with recent measurements from 
Song et al.l ll201ll); sligh tly lower than the estimates in 
Cabre fc Gaztafiagal | |200gh . 



on the scales as large as 200 h~ x Mpc. The effects of non- 
linear BAO diffusion become comparable to the error-bars 
on scales between 80 hT 1 Mpc and 200 h' 1 Mpc and should 
be taken into account to get accurate theoretical predictions. 
The effects due to a non-flat /i-distribution are next in order 
of importance after nonlinearities on large scales. For current 
data sets these effects are small compared to the statistical 
errors, but they will become important for future surveys. 
In what follows we will therefore only use the data on scales 
between 30 h' 1 Mpc and 200 hT 1 Mpc, and will ignore the 
nonlinear FOG, nonlinear growth other than BAO diffusion 
and wide-angle effects, but will take into account the effects 
of nonlinear BAO diffusion and the jj, distribution. 

In Sec. 13.71 we test the applicability of our model on 
the mock catalogs. The measurements of mean from the 
mocks have statistical error bars that are approximately nine 
times smaller compared to the SDSS data. To fit the mock 
measurements accurately we will also have to take into ac- 
count the FOG effect. 

To summarise, our theoretical model of the correlation 
function will be given by 



f* (r,fi,a)W(r,fi,a)P e (fi)dfida, 



(25) 



where the function £ th (r, fi, a) is computed by Fourier trans- 
forming a power-spectrum given by formula in Eq. (|10[) . We 
will model the real-space power-spectrum on the right hand 
side of Eq. (|10p as a linear power-spectrum damped with a 
Gaussian function of Eq. (|16[l to account for nonlinear dif- 
fusion of the BAO peak. The other effects are considered 
negligible for the SDSS data on these scales. 



3.7 Testing RSD models with mock catalogs 

To test our analysis of the effects that have to be taken into 
account to analyse RSD in SDSS DR7 data, and to esti- 
mate the statistical errors on our measurements (as shown 
on Fig. [2]), we use galaxy catalogs from the Large Suite of 
Dark Matter Simulations (LasDamas: ?jfl The LasDamas 
simulations are designed to model the clustering of Sloan 
Digital Sky Survey (SDSS) galaxies in a wide luminosity 
range and in the redshift range 0.16 < z < 0.44. The sim- 
ulations are produced by placing artificial galaxies inside 
dark matter halos using an HOD with parameters measured 
from the respective SDSS galaxy samples. We use 80 "Ori- 
ana" catalogs that have exactly the same angular mask as 
the SDSS survey and subsample them to match the redshift 
distribution of the Luminous Red Galaxies (LRG) in our 
SDSS DR7 data set. The LasDamas mocks have insufficient 
galaxies at redshifts below z < 0.2, as a result the mocks 
will slightly overestimate the shot noise. We do not expect 
this to be important since the effected region contains only 
a small fraction of the volume available. 

We apply exactly the same weighting to the mocks as 
to the real catalog and compute zeroth, second and fourth 
Legendre momenta of the redshift-space correlation function 
from them using Eq. J5). We also compute the normalised 
quadrupole Q(r) as given by Eq. (0. 

We estimate covariance matrices corresponding to the 



http : //lss .phy . vanderbilt . edu/lasdamas/ 
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statistical errors of our measurements, based on assuming 
that the Legendre momenta are drawn from a multi-variate 
Gaussian distribution 

cst at = 1_ J2 [X(r0 - X(n)] [X(r,-) - M^)] , (26) 

where X(r) is a vector of the measurements of at scale 
r for £ — 0, 2, 4 and X is the mean value from all 80 mock 
catalogs. 

The mean Legendre momenta measured from the Las- 
Damas mocks are shown in Fig. 1111 The error-bars corre- 
spond to the square root of the diagonal terms in the co- 
variance matrix C stat /80 and the lines show theoretical pre- 
dictions computed making different assumptions. Our theo- 
retical predictions, with the parameters of the simulations, 
provide a very good fit to the data. The bottom panel on 
Fig- HU shows that the theoretical prediction underestimates 
£4 on scales smaller than 50 /i -1 Mpc. The fourth Legendre 
moment measures a higher frequency fi dependence of cor- 
relation functions and therefore is more sensitive to different 
systematic effects. Since we are not using £4(7") in our fits we 
did not attempt to investigate this issue further. The mea- 
surements are strongly positively correlated and the error- 
bars presented here reflect only small part of the covariance 
matrix. 

Fiaure [T2l shows the measurements of Q from LasDamas 
mocks with the similar definition of the error-bars as Fig. 1111 
In the Kaiser formalism Q(r) is expected to be a straight 
line damped at smaller scales because of FOG effects. In 
reality the measured Q(r) will deviate from a straight line 
even within the Kaiser model at larger scales since it is es- 
timated by a discrete sum in Eq. instead of continuous 
integral. Figs.[ll]and[T2]clearly show that the model adopted 
in Sec. 13.61 can describe the measurements very well on all 
scales between 30 Mpc and 200 hT 1 Mpc, while using the 
Kaiser formula without modifications would fail to fit on 
scales around BAO peak and larger. 



4 TESTING COSMOLOGICAL MODELS 

Figure [2] shows that the signal to noise of the mea- 
sured £4 is very small on all scales so that it can 
not be used to extract RSD information. Consequently, 
for simplicity, we will not include measurements of £4 
in our analysis. Recent studies have shown that includ- 
ing hexadecapole in the fit improves errors on measured 
cosmological parameters IITaruva. Saito fc NishimichillioTll : 
iKazin. Sanchez &: Blantonll201ll ). We find similar improve- 
ment only if the geometrical information is measured from 
the survey itself through the AP effect without imposing 
any external priors. If the background geometry is fixed by 
strong external priors, as is the case in our analysis, the 
difference in measurements of bag and fog is very small. 

The normalised quadrupole Q by definition does not 
contain any extra information compared to £0 and £2- 
Figure [TU] shows that our measurements of Q are nois- 
ier than first two Legendre momenta on scales larger than 
50 h -1 Mpc. The analysis of Q measurements is in some way 
simpler, because the normalised quadrupole, under some 
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Figure 11. Measurements of mean |^| from 80 LasDamas "Ori- 
ana" mocks in a redshift range 0.16 < z < 0.44. Dashed line 
shows predictions of the Kaiser model (for the best fit values of 
berg, fag and <r v ), while solid line corresponds to the theoreti- 
cal predictions of the model with nonlinear BAO damping and 
non-flat /i-distribution (for the best fit values of berg , fcrg , a v and 
So). 
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Figure 12. Measurements of mean Q from 80 LasDamas "Ori- 
ana" mocks in a redshift range 0.16 < z < 0.44. Dotted line 
shows theoretical predictions of Kaiser model, dashed line shows 
predictions of the Kaiser model with FOG corrections, modeled 
with a Gaussian function, and accounting for the fact that Q is 
computed from a discrete sum, while solid line corresponds to the 
theoretical predictions of our model. Wiggles in solid and dashed 
lines are related to the fact that the integral in the definition of 
Q is replaced by discrete sum and we expect similar wiggles in 
the data. 



parameters A(z) describing the amplitude of the correlation 
function and its growth with redshift. Each model will also 
depend on the phenomenological parameter Eo describing 
nonlinear diffusion of the BAO peak. We will treat Eo as a 
nuisance parameter and marginalise over it with a uniform 
prior. For this reason, we do not include the Eo dependence 
of the likelihood in the equations given below. 

For the background expansion we will assume that the 
Universe is well described by a spatially-flat wCDM model 
composed of non-relativistic matter with energy density fi m 
some part of which is in baryons with energy density fib- 
The rest of the energy density, in this model, is assumed 
to be in a smooth dark fluid with the equation of state w. 
To complete the background model we need to specify the 
expansion rate of the Universe at present h — i/o/100, where 
Ho is a Hubble parameter. 

For the cosmological parameters describing the ob- 
served amplitude of clustering we will make three different 
assumptions ranging from the most specific model to more 
general assumptions. 

For every theoretical model we compute a x 2 function 

XL(P, A) = [X - X(r0]C tot [X - X(r J )] T , (27) 

where X(r) is a vector of the measured £o(r) and &(r), 
X(r) is the model to be tested, and the total covariance 
matrix is given by Eq. (|26|) . Assuming that the measurement 
errors are closed to Gaussian, the likelihood for a given set 
of cosmological parameters given data will be 



Ct 



exp(-Xtot/2). 



(28) 



general assumptions |f] does not depend on the shape of 
the power spectrum and only depends on the parameters 
describing the amplitude. The drawback is that with Q 
only a combination of growth and bias /3 = f/b can be 
measured, but not the two parameters individually, and 
Fig. 1121 sh ows that Q is model dependent (in detail see also, 
iTocchini-Valentini et afll201ll ). 

In our analysis we will perform a joint fit to £o and £2- 
We will explicitly use a prior based on the joint constraints 
from WMAP and SNIa data to deal with the ambiguity 
in the shape of the matter power-spectrum. This will al- 
low us to measure bias and growth at the same time and 
to extract more information from the correlation function. 
Fitting the first two even Legendre momenta of correlation 
function is advantageous compared with fitting to the two 
dimensional correlation function £(r) for two reasons: a pair 
of one-dimensional functions £0 and £2 are easier to visualise 
and work with and, as we showed above, they contain most 
of the cosmologically relevant information anyway; also the 
measurement errors on & are more Gaussian, as we show 
in Sec. 14.21 compared to the errors of f(r) and therefore the 
reconstruction of the likelihood surfaces is more robust. 

Our theoretical model of Legendre momenta of the cor- 
relation function will depend on a set of parameters p de- 
scribing background expansion of the Universe and a set of 



9 This does not hold, for example, when the fj, distribution does 
not correspond to an isotropic distribution of galaxy pairs. Non- 
linear effects will also make Q deviate from a straight line. 



4.1 Inaccuracies in the estimation of covariance 
matrix 

Estimating covariance matrices of galaxy two-point corre- 
lation function in configuration space is a nontrivial task. 
Many different techniques have been used before to tackle 
this issue, including internal procedures - b ased only on 
the observed da ta itself - such as jacknife jlLucevI ll979T ) 
and bootstrap (|Barrow. Bhavsar fc Sonoda Il984) meth- 
ods; analytical estimates of the errors ( Mo. Jing fc Bornerl 
1992); Monte-Carlo sampling of random initial conditions 
and combination of analy tical methods and Monte-Carlo 
(|Padmanabhan et al.ll2005T i. Studying the three-dimensiona l 
clustering on scales below 25 Mpc. lNorberg et al.l (|2009h 
showed that internal methods do recover the principal com- 
ponents of the real covariance matrix in a robust way but 
can not accurately reproduce the errors themselves, usually 
overestimating them by as much as 40 percent. 

Our statistical covariance matrices are estimated from 
the sample of 80 mock catalogs. This number is far less then 
sufficient to accura tely measure the errors and correlations. 
ICabre et al.l l|200ct ) used a sample of 1000 mocks, in their 
study of cross-correlation between the map of CMB tem- 
perature anisotropies and large scale structure, and found 
that even with 200 simulations the error bars could be un- 
derestimated by about 20 percent. 

One of the ways of reducing the effect of inaccurate co- 
variance matrix estimation is to find the eigenvectors of the 
normalized covariance matrix and then only use the eigen- 
modes that have high signal to noise, since there error esti- 
mates are expected to be more reliable. We do not attempt 
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to do this in our paper; this would remove large scale infor- 
mation which we are interested in. 

We will use the full covariance matrix estimated from 80 
mock catalogs as our best guess to the real structure of the 
measurement errors. This is good enough for the purposes 
of our current work. As we will show below the errors on 
current data are too big to result in tight constraints on 
cosmological parameters and we apply our method to the 
real SDSS-II data to provide a "proof of concept". Next 
generation of surveys, with significantly tighter error bars 
on the measurements of correlation function, will require a 
more thorough investigation of this issue. 

4.2 Inaccuracies in the posterior likelihood 
function 

Equation ([28} represents a true likelihood function only if 
the measurements of variables X that were used in comput- 
ing x 2 have errors that are distributed as multivariate Gaus- 
sian random variables. There are reasons to b elieve that the 
error s on fo, for instance, are not Gaussian (|Norberg et al.l 
2009). To check if the assumption of Gaussianity holds rea- 
sonably well for our measurements we take the measure- 
ments of £o and £2 at the different scales from all 80 mock 
LasDamas catalogs and construct normalized variables 

Y=^*, (29) 
0"x 

where X and ctx are the average value and dispersion com- 
puted from all 80 mocks. If the measured X are Gaussian, Y 
should be distributed according to the normal distribution 
with mean zero and variance one. 

Figs. 1 13l fc 1141 shows the distribution of Y for the mea- 
surements of £0 and £2- 

We perform a Kolmogor ov-Smirnov test (for deta iled 
discussion of the test see, e.g., ICorder fc Foreman! 1)20091) ) to 
see if the empirical distribution of measurements is con- 
sistent with the null hypothesis that they are drawn from 
a Gaussian distribution. The Kolmogorov-Smirnov test con- 
firms that the distribution of £0 is consistent with the null 
hypothesis at 13 percent confidence level and £2 is consistent 
at 25 percent confidence level. In both cases the small possi- 
ble deviations from Gaussian cumulative distribution func- 
tion reflect the fact that deviations above the mean value 
are slightly more likely then deviations below the mean at 
the tails of the distribution. 

To check how the Gaussianity of £e measurement er- 
rors depends on the scale we split r range into two with 
30h _1 Mpc < n < 75h _1 Mpc and 75h _1 Mpc < r 2 < 
200/i~ 1 Mpc and perform a similar Kolmogorov-Smirnov 
test on small scale and large scale measurements separately. 
Our empirical distribution of £0 is more consistent with the 
assumption of Gaussianity on small scales. For n the KS 
test accepts the null hypothesis at 36 percent, while for ri 
the null hypothesis is accepted at 8 percent. For the £2 the 
trend is opposite KS likelihood for n is 5 percent, while for 
r2 it is 17 percent. 

The variable log(l + £0) is a slightly better fit to the as- 
sumption of Gaussianity, with a KS likelihood of 23 percent 
over all scales. 

For our purposes the variables £0 and £2 are close 
enough to the Gaussian distributed variables and we con- 
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Figure 13. Histogram of the normalised scattering of go measure- 
ments and their cumulative distribution from 80 mock LasDamas 
catalogs compared to the normal distribution with mean zero and 
unit variance. 



elude that the usage of Eq. is justified for computing 
likelihood surfaces and confidence level intervals. 

The two dimensional redshift-space correlation func- 
tion £(<7,7r), where a and n are along the line-of-sight and 
across the line-of-sight separations, itself is often used in the 
analysis of RSD and BAO. We perform the same check on 
the measurements of this two dimensional correlation func- 
tion between the scales of 30-60 Mpc histogramed in 40 
two-dimensional bins. The resulting histogram is shown on 

Fig. [m 

Figure. \T5\ shows that the distribution of measured 
7r) around its mean value, at least for the mock cata- 
logs that we use, are not Gaussian. The null-hypothesis that 
these random variables are drawn from a Gaussian distribu- 
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Figure 14. Histogram of the normalised scattering of £2 measure- 
ments and their cumulative distribution from 80 mock LasDamas 
catalogs compared to the normal distribution with mean zero and 
unit variance. 



tion is rejected at more than 50 percent significance level. We 
conclude that if the \ 2 an d likelihood functions are denned 
as in Eqs. (|27|l and (I28|l . the Legendre momenta are more 
appropriate as variables compared to the two-dimensional 
correlation function. This is not surprising since the mea- 
sured Legendre momenta are the weighted sums of the two- 
dimensional correlation function and will always tend to be 
more Gaussian irrespective of the underlying distribution, 
thanks to the Central Limit Theorem. 
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Figure 15. Histogram of the normalised scattering of £(cr, it) 
measurements from 80 mock LasDamas catalogs. 

4.3 Covariance matrix as a function of 
cosmological parameters 

Another systematic effect in estimating covariance matrix 
is the dependence of C s at on the cosmological model. Our 
estimates of the covariance matrix are based on the mock 
catalogs that were created for a specific cosmological model, 
namely a spatially-flat ACDM with f2 m = 0.25, as = 0.8 
assuming that the gravity is described well by GR. In other 
cosmological models or different vales of parameters the in- 
trinsic scattering in the correlation function and therefore 
the covariance matrix will be different. The scaling of C stat 
with cosmological parameters is extremely difficult to model 
theoretically for nontrivial survey volumes. 

To estimate this effect we will again use the Fisher ma- 
trix calculations for an SDSS-II like survey. We compute a 
Fisher matrix F(b, f, an , a±_ , p) for different values of cosmo- 
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Figure 16. Fisher matrix predictions of the correlated errors 
of parameters b and fag from an SDSS-II like survey. Different 
contours correspond to the la confidence level ellipses for the 
fiducial cosmologies with fag equal to 0.2, 0.3, 0.4 and 0.5 from 
largest to smallest ellipse respectively. Bias is fixed to b = 2. 
Horizontal tickmarks shows the deviation from fiducial value. 



Figure 17. Fisher matrix predictions of the correlated errors 
of parameters bag and / from an SDSS-fl like survey. Different 
contours correspond to the la confidence level ellipses for the 
fiducial cosmologies with bag equal to 1.4, 1.5, 1.6 and 1.8 from 
outside to inside respectively. Growth rate is fixed to / = 0.45. 
Vertical tickmarks shows the deviation from fiducial value. 



logical parameters and look at how the expected errors on 
the measurements of growth scale. The uncertainties in the 
measurements of the power-spectrum can be schematically 
divided into two parts: coming from the cosmic variance and 
from the shot-noise. The shot-noise contribution depends 
on the total number of galaxies and their distribution in 
the survey volume and is insensitive to the underlying cos- 
mological model. The cosmic variance component depends 
on the parameters determining the overall amplitude of the 
power-spectrum b, f and ag, but is not very sensitive to the 
cosmological parameters describing its shape p. 

We derive Fisher matrix errors on the measurements of 
the growth parameter fag and bias bag for different fidu- 
cial values. This predictions are shown on Figs. 1161 and [P71 
Figure [16] shows the size of la ellipses for different fiducial 
values of fag when b = 2, while Fig. 1171 shows the same el- 
lipses for different values of bag when / = 0.45. The relative 
change is small compared to the sizes of the contours them- 
selves. We conclude that this effect is relatively unimportant 
for the range of values fag and bag allowed by our data and 
ignore it in our analysis. 

For next generation surveys, however, the errors on the 
measurements of growth will be significantly smaller and this 
effect will have to be taken into account. This implies that 
deviations from the best-fit value towards stronger clustering 
amplitude will be more likely then deviations of the same 
magnitude towards weaker clustering amplitude. 

4.4 wCDM and General Relativity 

In a specific cosmological model the growth rate will depend 
on the parameters describing background geometry as well 
as the theory of gravity. Assumptions about exact nature 



of this dependence bring in a very strong theoretical priors 
that might affect the results of data analysis strongly. In our 
work we will always make three separate assumptions about 
how the growth rate depends on the background expansion 
from most restricting to almost premise free. 

First we will assume that General Relativity (GR) is 
the correct theory of gravity. In this case the growth func- 
tion f ST (z) = f(p,z) can be computed at every redshift 
from basic cosmological parameters p. The other two num- 
bers that are necessary to completely describe the amplitude 
of the correlation function are the linear bias b(z) and the 
overall amplitude of clustering ag(z = 0). In this wCDM + 
GR model our cosmological parameters of interest will be 
p = (Q m ,h,w) and A = (bi(z),ag). 

We will use the prior likelihood on p and ag from the 
WMAP7 measurements and SNIa data. We do this by using 
the relaxed MonteCarlo Markov Chain of this joint datal 10 ! 
For every b(zi) we will go through the MCMC chain and for 
each value of p and ag compute the growth function f(z) 
and then the theoretical correlation function. Afterwards we 
will marginalize over p and ag. This is equivalent to taking 
the following analytical integral 

C £I (b(z l )) = J £ t ot(p, A)£ prior (p, a s )dpdag, (30) 

where £ pr ior is effectively given by the MCMC. This will 
enable us to derive constraints on the linear bias parameter 
b(z) in two redshift bins. 
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4.5 7 parametrization of growth 

For our sec ond model w e will consider the 7 parametrization 
of growth (|Linderll2005l V In this model the growth function 
is assumed to depend on parameters p as 



/(*) = 



041 + z) 3 



0»(1 + zf + (1 - Om)(l + z)~^+™) 



(31) 



where 7 is a redshift and scale independent number. For 
7 = 0.55, Eq. (|31[) gives numerical results that are very 
close to the predictions of GR. If 7 is larger than 0.55 the 
growth is weaker compared to GR and vice versa. 

In this model A = (b(Zi), 7, as). We will assume that 
the shape of the correlation function can still be accurately 
modelled by wCDM predictions and will use WMAP + SNIa 
MCMC to marginalize over p and as as in Sec. 14.41 

To use the WMAP prior on as for 7 parametrization we 
have to take into the account the fact that erg is not directly 
measured by CMB experiments. The measured quantity is 
as at the last scattering surface with z ~ 1100 and then the 
as(z = 0) is inferred by rescaling 



a 8 (z = 0) = a 8 (z = 1100) 



G(z = 0) 
G(z = 1100) ' 



(32) 



where G(z) is the growth factor assuming GR. 

To make the priors on as as given by WMAP MCMC 
chain consistent with our assumption that the growth of 
structure is modified as in Eq. (|31[) we rescale the values of 
as in WMAP MCMC chain by 



: G(z = 1100) GR G(z = 0) 



a 8 (z = Of = a 8 (z = Of 



(33) 



G(z = 1100)t G(z = 0) gr: 

where G 7 and G GR are the growth functions computed in 7 
parametrization and GR respectively. After marginalization 
we will get a posterior likelihood function C 1 {b(zi) , 7) . The 
measurements of 7 in general will be correlated with the 
measurements of bias. 



4.6 Free growth 

In the last model we will not make any assumptions about 
the relationship between / and p and will treat /(z;) in each 
redshift bin as a free parameter. In this model three parame- 
ters describing the amplitude /, b and as are degenerate and 
only two combinations of them can be measured indepen- 
dently. We will choose these combinations to be b(zi)as(zi) 
and f(z l )a a (z i ). 

We will assume again that only the growth of the per- 
turbations is different from the GR case and the overall 
shape of the correlation function can still be model by 
wCDM model. If we make this assumption we can use the 
same chains to marginalize over p so that we are left with the 
posterior likelihood function £ ls (b(zi)as(zi), f{zi)as{zi)). 



5 RESULTS AND DISCUSSION 



We use the method outlined in Sec. [3] to constrain pa- 
rameters describing the redshift evolution of the clustering 
of LRGs. We first use only scales up to 60 h~ x Mpc and then 



use all scales up to 200 h 1 Mpc. Our results are presented 
in Tab. Q] 

For the wCDM + GR model the small scale data con- 
strains the real-space amplitude of the galaxy clustering 
signal in both redshift bins with the accuracy of about 5 
percent. These measurements are consistent with previous 
estimates, showing that LRGs are highly biased tracers of 
the underlying matter field. The power spectrum amplitudes 
in the two redshift bins are close and consistent with the 
assumption of the constant clustering amplitude. The con- 
straints improve when we extend the analysis to include data 
on scales 60 - 200 /i -1 Mpc. This inclusion results in slightly 
higher estimates of bias in higher redshift bin, but the two 
measurements are consistent at a la confidence level. 

For the more general 7 parametrization that has GR 
as a specific case, the 7 parameter is constrained with the 
precision of about 22 percent with small scale data and 19 
percent when larger scales are included. Both large and small 
scale measurements prefer a weaker growth than in GR but 
are consistent with GR results at la confidence level. When 
the large scale data is included the best-fit values for 7 are 
closer to the GR values. The 7 parametrization fits give best- 
fit values of bias in both redshift bins that are consistent 
with those measured assuming GR, following the standard 
degeneracy between the bias and the RSD signal. 

Note that when adopting the 7 parametrization, we 
have implicitly assumed that the growth is modified with 
respect to GR in a manner that is independent of the scale. 
Even if the real modifications of gravity are scale dependent, 
the 7 parametrization will still be able to capture deviations 
from GR, but the measured 7 will be an average over the 
scales being considered. 

For the most general model of free growth the param- 
eter fas can be measured with the accuracy of about 15 
percent in both redshift bins. The inclusion of large scale 
data, again improves these constraints slightly. This shows 
that at larger scales SDSS DR7 clustering data is noisier 
and introduces more scatter. The best fit values of growth, 
when it is allowed to freely vary are consistent with the pre- 
dictions of GR. Recovered values of bias are consistent with 
the ones measured in previous two models. These measure- 
ments should be considered as giving the average growth 
over the scales covered. 

Similar measurements for bas and fas based on fits 
to measured two-dim ensional £ s (r,n) have been made be- 
fore fr om SDSS DR6 (ICabre fe Gaztanagall2009l ) and SDSS 
DR7 (|Song et al.llioilF based catalogs. Our constraints on 
bias (about 5 percent in both redshift bins) are comparable 
to previously derived results (about 6 percent at lower red- 



shift a nd 8 percent at higher redshift in Cabrc & Gaztanaga 



(2009); about 3 percent in both redshift ranges in Song et al 



(|201lj)), while our constraints on fas (about 16 percent 
at lower redshift and 15 percent at higher redshift) are 
stronger than results derived in previous studies (19 per- 
ce nt at lower redshifts and 2 2 percent at higher redshifts 
in ICabre fc Gaztanagal ([2009) ; 15 perc ent at lower reds h- 
fits and 12 percent at higher redshift in ISong et al.l (|201ll )L 
Slightly improved sensitivity to fas could be due to the fact 
that we are fitting to Legendre momenta function rather 
than the two-dimensional correlation function itself. 

The measurements of the pair of variables bas , fas and 
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Model 


Variable 


Scales less than 60 Mpc//i 


Scales up to 200 Mpc/h 


"Standard" model expectation 


wCDM 


b(zi)as(zi) 
£>( 2 2)°-8(>2) 


1.4216 ± 0.0724 
1.4053 ± 0.0582 


1.3890 ± 0.0448 
1.50565 ± 0.0352 




7 


b{zx)cr%(z\) 

7 


1.4641 ± 0.0790 
1.4641 ± 0.0703 
0.7366 ± 0.1638 


1.4156 ± 0.0491 
1.5107 ± 0.0395 
0.5842 ± 0.1116 


0.55 


Free growth 


6(zi)ct 8 (2 1 ) 
b(z2)cr$(z2) 

/(Zl)<7 8 (zi) 
/(2 2 )o- 8 (22) 


1.4663 ± 0.0828 
1.4511 ± 0.0758 
0.3665 ± 0.0601 
0.4031 ± 0.0586 


1.4157 ± 0.0521 
1.5092 ± 0.0398 
0.3512 ± 0.0583 
0.4602 ± 0.0378 


0.4260 
0.4367 



Table 1 . Constraints on parameters describing growth and clustering bias of galaxies with respect to the matter field in different models 
with and without including measurements from scales more than 60 /i -1 Mpc. "Standard" model refers to the spatially-flat ACDM with 
f^m, eg = 0.8 and general relativity. 



6<rg, 7 are very weakly correlated and can be assumed to be 
independent for all practical purposes. 



6 CONCLUSIONS 

In this paper we have considered systematic deviations 
from the linear plane-parallel RSD model for the large 
scale clustering of galaxies. These include systematic devi- 
ations due to wide-angle and non-linear effects, and prob- 
lems caused by inaccurate modelling of the redshift dis- 
tribution. By testing different models against the measure- 
ments from N-body simulations we checked that, by includ- 
ing these effects, we can fit simulated large-scale RSD data 
extremely well. We have also considered the relative impor- 
tance of these effects, showing that the wide-angle effects 
are small for the SDSS DR7 survey and can be safely ig- 
nored even on scales as large as 200 h' 1 Mpc, but nonlinear 
damping of the Baryon Acoustic Oscillation (BAO) peak 
l|Meikson. White fc Peacock! 1 19991 ) has to be taken into ac- 
count in order to properly fit the data. The non-isotropic /1- 
distribution effects are small for SDSS DR7 but much larger 
compared to wide-angle effects and will be important for 
future surveys. 

Currently available measurements of galaxy clustering 
are low signal to noise on very large scales and the inclu- 
sion of systematics that we discussed above do not bias the 
estimates of growth and bias. For the next generation of 
surveys (e.g., BOSS, EUCLID) the effects of uncertainties 
in radial selection, large-scale nonlinearities and non-flat \i- 
distribution will be comparable or even larger than statisti- 
cal errors. 

We did not account for FOG effect in our fits. For 
the spherically averaged correlation function FOG is ex- 
pected to be small on the scales that we consider but for 
the quadrupole of correlation function this effect could be 
comparable to the errorbars on the scales up to 40 K~ Y Mpc 
(See, Fig. [9} Since we use Eq. (|16|l . which already includes 
FOG effects, to model large scale nonlinearities, applying ad- 
ditional FOG damping term would not be consistent. For the 
analysis of future high precision measurements of clustering, 

11 For Q the effect is very large even on large scales since the 
normilized quadrupole depends on the integral over correlation 
function over all scales. 



proper modelling of FOG effects on small and intermediate 
scales will be necessary. 

In our analysis we only kept linear order terms in real- 
to-redshift space mapping. Recent works have demonstrated 
that nonlinear contribution to Eq. (JH| introduce additional 
bias in theo retical estimates of Pes and Pse and there- 
fore in &(r) dTaruya. Saito fc Nishimichill201ll ; lTandl201ll ; 
iReid fc White] 201 ll ). Future measurements will also require 
careful treatment of these nonlinear effects. 

Different ways of extracting RSD information have been 
considered before, including a fit to the two-dimensional cor- 
relation function £(cr, n) and the normalised quadrupole Q. 
We argue that the best approach is to perform a joint fit 
to measured Legendre momenta of the correlation function. 
Based on the simple linear plane-parallel model they contain 
exactly the same information as £(cr, n) and their measure- 
ment errors are more Gaussian, which makes the interpreta- 
tion of straightforward. Compared to only using the nor- 
malised quadrupole Q, they contain significantly more infor- 
mation and allow for the measurements of bias and growth 
independently instead of measuring only their ratio. Also 
Q was originally proposed because it was believed to have 
certain advantages of being independent of the shape of the 
power-spectrum and nonlinear effects. We show that that 
does not hold on large scales: Q is affected by nonlinearities 
as much as the correlation itself and it is still affected by 
AP effect and the dependence on the shape of the power- 
spectrum and on the background cosmological parameters 
is not completely removed. 

We have analysed the SDSS DR7 LRG clustering in 
redshift-space and obtained constraints on bias and param- 
eters describing structure growth in two redshift bins. We 
have presented what we consider to be a very robust analy- 
sis, taking into account all of the effects that could influence 
the redshift-space correlations function. The inclusion of the 
very large scale data does not improve our measurements of 
bias and growth parameters significantly: current measure- 
ments of the correlation function on the scales larger than 
60 h' 1 Mpc are too noisy to be of practical interest. The next 
generation of ongoing and planned surveys such as BOSS; 
ISchlegel, White fc Eisensteinl 120091') BigBOSS and the ESA 
Euclid mission ( Laureiis et a.l.ll2009T ) will enable us to mea- 
sure clustering properties of galaxies on very large scales 
with high accuracy. For these surveys where the measure- 
ments are more precise, the full treatment of RSD effects 
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will be very important. In addition, there will be signifi- 
cantly more information available on the largest scales. 



ACKNOWLEDGEMENTS 

We thank anonymous referee whose comments and sugges- 
tions helped us to improve our paper. LS is grateful for sup- 
port from the European Research Council. WJP is grateful 
for support from the UK Science and Technology Facilities 
Council, the Leverhulme trust and the European Research 
Council. AR is grateful for the support from a UK Science 
and Technology Facilities Research Council (STFC) PhD 
studentship. We thank Beth Reid for useful comments on 
an early draft of this paper and Marc Manera for useful 
discussions. LS acknowledges partial support from Georgian 
National Science Foundation grant GNSF ST08/4-442 and 
SNFS SCOPES grant no. 128040. We are grateful to Las- 
Damas project for making their mock catalogs publicly avail- 
able. Numerical computations were done on the Sciama High 
Performance Compute (HPC) cluster which is supported by 
the ICG, SEPNet and the University of Portsmouth and 
COSMOS consortium supercomputer within the DiRAC fa- 
cility jointly funded by STFC, the Large Facilities Capital 
Fund of BIS and the University of Portsmouth. 



REFERENCES 

Abazajian, K. N., et al., 2009, ApJS, 182, 543 
Alcock C, Paczynski B., 1979, Nature, 281, 358 
Ballinger W. E., Peacock J. A., Heavens A. F., 1996, MN- 
RAS, 282, 877 

Barrow J. D., Bhasvar S. P., Sonoda D. H., 1984, MNRAS, 
210, 19 

Bharadwaj S., 1996, ApJ, 472, 1 

Cabre A., Fosalba P., Gaztanaga E., Manera M., 2008, MN- 
RAS, 381, 1347 
Cabre A., Gaztanaga E., 2009, MNRAS, 393, 1183 
Carlson J., White M., Padmanabhan N., 2009, PRD, 80, 
043531 

Cole S., Fisher K.B., Weinberg D.H., 1995, MNRAS, 275, 
515 

Corder G. W., Foreman D. I., 2009, Nonparametric statis- 
tics for non-statisticians: A step-by-step approach, Wiley- 
Blackwell, Hoboken, NJ 
Crocce M., Scoccimarro R., 2006, PRD 7 3, 063519 
Colless M., et al., 2003, preprint astro-ph/0306581 
Crocce M., Scoccimarro R., 2008, PRD77, 023533 
Crocce M., Gaztanaga E., Cabre A., Carnero A., Sanchez 
E., 2011, MNRAS, 417, 2577 
Dalai N., Dore O., Huterer D., Shirokov A., 2008, PRD, 
77, 123514 

Desjacques V., Seljak U., 2010, Adv. Astron., 2010, 908640 
Desjacques V., Seljak U., Iliev IT., 2010, MNRAS, 396, 85 
Eisenstein D.J., et al., 2001, AJ, 122, 2267 
Eisenstein D.J., Seo H.-J., White M., 2007, ApJ, 665, 14 
Feldman H.A., Kaiser N., Peacock J.A., 1994, ApJ, 426, 23 
Feldman H. A., Watkins R., Hudson M. J., 2010, MNRAS, 

407, 2328 
Fisher K.B., 1995, ApJ, 448, 494 



Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shi- 
masaku, K., Schneider, D. P., 1996, AJ, 111, 1748 
Gunn, J. E., et al., 1998, AJ, 116, 3040 
Guzzo L. et al., 2009, Nat., 541, 451 

Hamilton A.J.S., 1992, ApJ, 385, L5 

Hamilton A.J.S., 1997, preprint ( |arXiv:astro-ph/9708102j | 
Hawkins E., et al, 2003, MNRAS, 346, 78 
Jackson J.C., 1972, MNRAS, 156, 1 

Jennings E., Baugh C. M., Pascoli C, 2011, ApJ, 727, L9 
Jing Y.P., Borner O, 2004, ApJ, 617, 782 
Kaiser N., 1987, MNRAS, 227, 1 

Kashlinsky A., Atrio-Barandela F., Kocevski D., Ebeling 

H., 2009, ApJ, 686, L49 
Kazin E., et al., 2010, ApJ, 710, 1444 

Kazin E., Sanchez A.G., Blanton M.R., 2011, MNRAS, in 
press 

Landy, S. D., Szalay, A. S., 1993, ApJ, 412, 64 
Laureijs R. et al., 2009, Euclid Assessment Study Re- 
port for the ESA Cosmic Visions, (arXiv:0912.0914 [astro- 
ph.CO]) 

Lewis A., Challinor A., Lasenby A., 2000, ApJ, 538, 473 
Li C, Jing Y.P., Kauffmann O, Borner G., Xi K., Wang 

L., 2007, MNRAS, 376, 984 
Linder E., 2005, PRD, 72, 043529 
Lucey J. R., 1979, PhD thesis, Univ. Sussex 
Macaulay E., Feldman H., Ferreira P. G., Hudson M. J., 

Watkings R., 2011, MNRAS, 414, 621 
Matsubara T., 2004, ApJ, 615, 573 
Matsubara T., 2008a, PRD, 77, 063530 
Matsubara T., 2008b, PRD, 78, 083519 
Meikson A., White M., Peacock J. A., 1999, MNRAS, 304, 

851 

Mo H. J., Jing Y. P., Borner G., 1992, ApJ, 392, 452 
Norberg P., Baugh CM., Gaztanaga E., Croton D.J., 2009, 

MNRAS, 396, 19 
Nusser A., Davis M., 2011, ApJ, 736, 93 
Okamura T., Jing Y. P., 2011, ApJ, 726, 11 
Okumura T., Matsubara T., Eisenstein D.J., Kayo I., Hik- 

age C, Szalay A.S., Schneider D.P., 2008, ApJ, 676, 889 
Padmanabhan N., et al., 2005, PRD, 72, 043525 
Papai P., Szapudi I., 2008, MNRAS, 389, 292 
Peacock J. A., Dodds S.J., 1996, MNRAS, 280, L19 
Percival W. J. et al., 2004, MNRAS, 353, 1201 
Percival W. J. et al., 2010, MNRAS, 401, 2148 
Press W.H., Teukolsky S,A., Vetterling W.T., Flannery 

B.P., 1992, Numerical recipes in C. The art of scientific 

computing, Second edition, Cambridge: University Press. 
Raccanelli A., Samushia L., Percival W. J., 2010, MNRAS, 

409, 1525 

Reid B. A., et al., 2010, MNRAS, 404, 60 

Reid B. A., White M., 2011, MNRAS, 417, 1913 

Ross A. J., Percival W. J., Crocce M., Cabre A., Gazganaga 

E., 2011, MNRAS, 415, 2193 
Samushia L., et al. 2010, MNRAS, 410, 1993 
Sanchez A., Crocce M., Cabre A., Baugh CM., Gaztanaga 

E., 2009, MNRAS, 400, 1643 
Schlegel D., White M., Eisenstein D.J., 2009, The astron- 
omy and astrophysics decadal survey, science white paper 
No. 314, ( arXiv:0902.4680l [astro-ph.CO]1 
Schlegel et al., 2009, preprint l|arXiv:0904.0468l [astro- 
ph.CO]) 

Scoccimarro R., 2004, PRD, 70, 083007 



© 0000 RAS, MNRAS 000, 000-000 



20 L. Samushia et al. 



Seljak U., Hamaus N., Desjacques V., 2009, PRL 103, 
091303 

Simpson F., Peacock J., 2010, PRD, 81, 043512 

Smith R. E., Peacock J. A., Jenkins A., White S. D. M., et 

al., 2003, MNRAS, 341, 1311 
Song Y.-S., Sabiu C, Kayo I., Nichol B., 2011, JCAP, 05, 

020 

Szalay A. S., Matsubara T., Landay S. D., 1998, ApJ, 498, 
LI 

Szapudi I., PRD, 70, 083536 
Takada M., PRD, 2006, 74, 043505 
Tang J., Kayo I., Takada M., 2001, MNRAS, 416, 229 
Taruya A., Nishimichi T., Saito S., 2010, PRD, 82, 063522 
Taruya A., Saito S., Nishimichi T., 2011, PRD, 83, 103527 
Tegmark M. et al., 2006, PRD, 74, 123507 
Thomas S. A., Abdalla F. B., Lahav O., 2011, PRL, 106, 
241301 

Tocchini-Valentini D., Barnard M., Bennett C. L., Szalay 
A. S., 2011, preprint l|arXiv:1101. 26081 [astro-ph.CO]) 

Watkins R., Feldman H. A., Hudson M. J., 2009, MNRAS, 
392, 743 

York D.G., et al., 2000, AJ, 120, 1579 

Zaroubi S., Hoffman Y., 1993, preprint 

( arXiv:astro-ph/9311013| 
Zehavi I. et al., 2005, ApJ, 621, 22 



© 0000 RAS, MNRAS 000, 000-000 



